Convection in Binary Fluid Mixtures. 
I. Extended Traveling Wave and Stationary States. 
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Nonlinear, spatially extended structures of convection rolls in horizontal layers of binary fluids 
heated from below are investigated in quantitative detail as a function of Rayleigh number for several 
negative and positive Soret coupling strengths (separation ratios) and different Lewis and Prandtl 
numbers characterizing different mixtures. A finite difference method was used to solve the full 
hydrodynamic field equations in a vertical cross section perpendicular to the roll axes, subject to 
realistic horizontal and laterally periodic boundary conditions in a range of experimentally accessible 
parameters. We elucidate the important role that the concentration field plays in the structural 
dynamics of the nonlinear states of stationary overturning convection (SOC) and of traveling wave 
(TW) convection investigated here. Structural differences in the concentration boundary layers and 
of the concentration plumes in TW's and SOC's and their physical consequences are discussed. These 
properties show that the states considered here are indeed strongly nonlinear, as expected from the 
magnitude of advection and diffusion in the concentration balance. The bifurcation behaviour of the 
states is analysed using different order parameters such as flow intensity, Nusselt number, a newly 
defined mixing parameter characterized by the variance of the concentration field, and the TW 
frequency. For further comparison with experiments, light intensity distributions are determined 
that can be observed in side-view shadowgraphs done with horizontal light along the roll axes. 
Furthermore, detailed structural analyses of all fields are made using colour coded isoplots, vertical 
and lateral field profiles, and lateral Fourier decompositions. They reveal, among other things, that 
the mirror-glide operation of lateral translation by half a wavelength combined with vertical reflection 
at the horizontal midplane of the layer is the longest persistent symmetry of TW and SOC states. 
Transport properties of the TW state are also discussed, in particular the mean lateral concentration 
current that is caused by the phase difference between concentration wave and velocity wave and 
that is roughly proportional to the TW frequency. This current plays an important role in the 
structural dynamics and stability of the spatially-localized traveling-wave convection investigated in 
an accompanying paper. 
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I. INTRODUCTION 

Thermally-driven convection in a binary fluid of two 
miscible components like ethanol- water or 3 He- 4 He has 
attracted much research activity lately. It has all the 
qualities of a paradigmatic system for studying prob- 
lems related to hydrodynamic instabilities, bifurcation 
theory, structure formation, complex spatio-temporal be- 
havior, and turbulence [1]: The experimental realization 
is sufficiently simple to allow for controlled experiments, 
and the governing field equations are well known. Fur- 
thermore, experiments [2-45] and theoretical approaches 
[46-77] that capture properties of the experimental set- 
ups, albeit to quite a varying degree, have observed an 
enormous richness of spatio-temporal phenomena. Here 
we investigate, by numerical integration of the hydrody- 
namic field equations, spatially-periodic, extended con- 
vective structures consisting of straight rolls. They ap- 
pear either as stationary patterns that we call states of 
stationary overturning convection (SOC) rolls for histor- 
ical reasons [2] or as traveling-wave (TW) patterns of 
propagating rolls. The analysis of these extended states 
is a useful and indeed necessary basis for an understand- 
ing of more complex structures like localized traveling 
wave (LTW) states, which we investigate in an accompa- 
nying paper [46]. 

In particular, we want to elucidate the role of the con- 
centration field. It is the combination (i) of the linear 
coupling of its degrees of freedom to the temperature 
field via the Soret effect and (ii) of its contribution to 
the buoyancy force — the concentration field directly in- 
fluences the driving mechanism for convection — that 
causes the additional richness of convective dynamics 
and structures in mixtures in comparison to pure one- 
component fluids. The latter exhibit as a first instability 
just a forwards bifurcation of the quiescent conductive 
state into stationary rolls. On the one hand, concentra- 
tion is advected, and since the concentration diffusion 
constant is small it is transported almost passively by 
the flow, except for boundary layer effects. But on the 
other hand the structural dynamics of the concentration 
field actively feeds back via the associated density varia- 
tions into the buoyancy force that drives convection. The 
Soret coupling constant -0 can be positive or negative de- 
pending on mean temperature and concentration, and 
for room temperature ethanol-water mixtures ^-values 
between about —0.5 and +0.2 can be easily realized ex- 
perimentally [3]. Without Soret coupling, -0 = 0, con- 
centration deviations from the mean eventually diffuse 
away. For 0^0, however, the externally-imposed tem- 
perature stress sustains concentration variations against 
the mixing action of convection via the Soret effect; these 
in turn alter the driving force for convection. For neg- 
ative Soret coupling, ip < 0, this interplay causes the 
appearence of a Hopf bifurcation into a pair of symmetry 
degenerated left and right TW states in addition to the 
bifurcation of SOC rolls. Both of these convective so- 



lution branches typically show a subcritical bifurcation 
topology — the bifurcation is forwards only in the im- 
mediate vicinity of the pure-fluid limit of vanishing Soret 
coupling, -0 = 0. The associated discontinuous onset of 
convection reflects the nonlinear feed back mechanism: 
convective homogeneization of the concentration distri- 
bution reduces the adverse buoyancy-depressing Soret- 
induced conductive concentration gradient and thus am- 
plifies buoyancy and convection intensity. 

We have calculated and analysed various states for 
negative as well as for positive Soret coupling param- 
eters. However, our main interest has been focussed 
on the experimentally observable TW and SOC states 
that are located at ip < on the upper solution branch 
of the bifurcation curves of, say, flow amplitude versus 
Rayleigh number. We should like to stress that these 
experimentally observable convective states are strongly 
nonlinear. Their spatio-temporal behaviour can neither 
be described quantitatively by a weakly nonlinear ex- 
pansion around their bifurcation thresholds nor by sim- 
ple heuristic amplitude equation approaches. We have 
therefore solved the full hydrodynamic field equations 
numerically with realistic boundary conditions using an 
explicit finite-difference method. We restricted ourselves 
to 2d convection in the form of straight rolls by suppress- 
ing spatial variations along the roll axes. In most cases, 
we enforced the wavelength of the states by imposing lat- 
erally periodic boundary conditions with the cell length 
equal to twice the height of the fluid layer, i.e., close to 
the wavelength seen in experiments. 

Since a review [47] of the early literature on convection 
in binary mixtures and a critical appraisal [1] of many of 
the more recent works is available - see also [48,49] - we 
add here only some short complementary comments that 
are mostly related to experimental and theoretical devel- 
opments relevant to our work. After the first unambigu- 
ous visual observation [4] of TW's in an ethanol-water 
mixture, the use of annular convection channels [5,6] 
that (i) eliminated endwall effects and (ii) enforced lat- 
eral periodicity facilitated comparisons with theoretical 
and numerical approaches. However, on the theoretical 
side, amplitude-equation models and few-mode Galerkin 
approximations (even after strongly increasing the num- 
ber of modes [7]) have failed to properly describe nonlin- 
ear TW's, e.g., their frequency variation with Rayleigh 
number, their bifurcation behaviour on the upper stable 
branch, and their field structure. So far, only Bensimon 
et al. [50-52] succeeded in calculating TW's on the up- 
per branch with a saddle connection to the lower back- 
wards bifurcating branch using a boundary layer model 
[78] for the concentration field. But this approach also 
poses some problems. 

On the numerical side, simulations of so called thermo- 
solutal convection by Deane et al. [53] yielded TW and 
SOC solutions. However, they are not fully applicable 
for quantitative comparison with the experiments [2-45] 
in binary mixtures for several reasons [1]. Stable non- 
linear TW's comparable to the experimental ones [5,6] 



have been obtained numerically in Ref. [54] using a fi- 
nite difference method. This analysis [54] showed the 
existence of large-scale mean Eulerian lateral currents of 
heat and concentration and a nonzero mean velocity field 
that had previously been found in few-mode Galerkin 
models by Linz et al. [55]. Also, the Lagrangian mo- 
tion in passive tracer transport [8,9] in TW states is now 
well understood [55,56,79,57,58,10] and is different from 
the above mentioned Eulerian currents and fields. Us- 
ing a finite-difference method at ip = —0.11, Yahata [59] 
calculated an SOC state on the upper solution branch 
showing boundary layer behaviour in the concentration 
field. A comparison with a 52-mode Galerkin model re- 
vealed that temperature and velocity fields could very 
well be described by few modes, but that the concen- 
tration field was not reproduced properly by this mode 
ansatz. Later, Yahata [60] calculated oscillatory convec- 
tion that was spatially confined to the vicinity of a lat- 
eral side wall. He also presented the concentration field, 
which had a similar structure as in an infinite system. 

Analysing experiments in an annular channel at ip = 
—0.25, Ohlsen et al. [11] found good agreement with 
the concentration boundary layer model of Bensimon 
et al. [51] and with numerical simulations [57] of the 
TW states on the upper branch all the way between the 
saddle and the nonhysteretic transition TW — > SOC. In 
short, straight rectangular channels, the lateral bound- 
aries make this transition hysteretic with a jump in TW 
frequency [12,13]. A comparison of experimental top- 
view shadowgraph images of TW and SOC states [14] 
with numerical results [57] also allowed the identification 
of the contribution from the concentration field to the 
light intensity profile caused by variations in the refrac- 
tive index of the mixture. This seems to have been the 
first time that structural properties of the concentration 
field - albeit vertically averaged - in binary fluid convec- 
tion have been observed experimentally. Subsequently, 
Winkler and Kolodner [15] extracted more direct infor- 
mation on the lateral structure of the concentration dis- 
tribution in TW and SOC states from side-view shadow- 
graph images obtaining good agreement with numerical 
results [57]. On the other hand, the shadowgraph mea- 
surements of Liu and de Bruyn [16], done with very nar- 
row (0.3d in ^-direction) and short rectangular channels, 
showed significant quantitative differences in intensities 
and bifurcation behavior with experiments in broader an- 
nular channels [14,15] and with our 2d numerical cal- 
culations. The differences can be traced back in large 
part to the cell geometry. In fact, our calculations agree 
well with experiments in channels of moderate width (1.3 
to 3 times the height) [15,14]. In very broad channels, 
straight rolls undergo 3d instabilities [80,61,17,18]. In 
very narrow ones, the influence of the walls becomes too 
strong [16]; in addition, the critical temperature differ- 
ence between the horizontal plates increases so strongly 
that non Oberbeck-Boussinesq effects become important. 
We should mention also that Zimmermann et al. [19-21] 
did side-view differential interferometry in ethanol-water 



mixtures. They found that, for mean fluid temperatures 
of about 10°C, the existence range of stable TWs is 
strongly reduced and at a mean temperature of about 
5°C, there were no stable TWs. 

Extended convective states also have been investigated 
for positive Soret coupling ip. For a survey, see, e.g., 
[22,23,62]. Lhost and Platten [24] observed, in a nar- 
row rectangular cell close to onset, very big rolls, as one 
would expect from a linear stability analysis which yields 
a vanishing critical wavenumber. Various experiments 
[25,22,26,13,23] done at not-too-small positive tp in broad 
rectangular and cylindrical cells found stationary, square- 
shaped convective structures. In agreement with the 
amplitude-equation model of Clune and Knobloch [62], 
they occur in the so-called Soret regime between the on- 
set r s tat and the reduced Rayleigh number r = 1 at which 
SOC rolls bifurcate in pure fluids. However, the channel 
walls of narrow straight or annular cells suppress square 
patterns and favour SOC rolls with axes oriented per- 
pendicular to the walls. In the so-called Rayleigh region 
above r « 1, square convection loses stability to a sta- 
ble oscillatory state [26,63,23]. Miiller and Liicke [63] ob- 
tained such temporal oscillations between rolls of perpen- 
dicular orientation, using a 10-mode Lorenz model with a 
phenomenological sidewall forcing. With further increase 
of the Rayleigh number, convection in the mixture ap- 
proaches the behavior in a pure fluid, i.e., SOC rolls of 
one direction are stabilized. To study square convection 
and the above described oscillatory state, Bestehorn [64] 
has numerically integrated the 3d field equations in the 
limit of infinite Prandtl number, albeit with very mod- 
erate spatial resolution of the concentration boundary 
layers. 

Our work is organized as follows: Section II describes 
the system we study, the field equations, the role of the 
Soret coupling, and the bifurcation thresholds of con- 
vective solutions from the quiescent conductive state. In 
section III, we analyse, for negative Soret coupling ip, the 
bifurcation behaviour of nonlinear TW and SOC states 
and their field structure and symmetry properties, espe- 
cially focussing on the concentration field and the trans- 
port properties. In section IV, we discuss how structure 
and bifurcation behaviour change with ip, and section V 
presents stationary states for positive Soret coupling. In 
section VI, we investigate variations with changing Lewis 
and Prandtl numbers. The conclusion in section VII lists 
our main results on nonlinear extended SOC and TW 
states. In the accompanying paper [46], we investigate 
spatially-localized traveling-wave convection. 



II. THE SYSTEM 

We consider a horizontal layer of a binary fluid mix- 
ture like alcohol- water or normal 3 He- 4 He under a ho- 
mogeneous gravitational field, g = — g e z , that is di- 
rected downwards. A positive temperature difference, 
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-upper ' 



between the lower and upper con- 



fining boundaries is imposed externally, e.g., via highly- 
conducting plates in experiments. When the resulting 
buoyancy force (cf. below) exceeds a threshold, convec- 
tive motion sets in. We consider here convection in the 
form of straight parallel rolls as seen in many experi- 
ments. Ignoring variations along the roll axes, we de- 
scribe 2d convection in an x — z plane perpendicular to 
the axes. 



A. Equations 

The system is described by the balance equations for 
mass, heat, concentration, and momentum in Oberbeck- 
Boussinesq approximation [81-83,47] 
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d t SC = -V • J ; 

J = uSC - LV (SC -tpST) , 

d t u = -V (u:u+p-crV:u)+B 

B = aR(ST + SC)e z . 
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Lengths are scaled with the height d of the layer, the time 
with the vertical thermal diffusion time cP/k, and the ve- 
locity field u = (u, 0, w) with n/d. Here, k is the thermal 
diffusivity of the mixture. ST = (T - T ) / AT denotes 
the deviation of the temperature from the mean temper- 
ature Tq in the fluid scaled by the temperature difference 
between the plates. Some of the unsealed quantities are 
underlined for distinction from the scaled ones. The field 
SC = (C - Cjo) (3 1 (a A T) is the scaled deviation of the 
mass concentration C_ = p 1 /(p 1 + p 2 ) of the solute from 
its mean Cq. Here p 1 and p 2 are the mass density fields of 
the two components. The thermal expansion coefficient 
a and the solutal expansion coefficient (3 of the fluid come 
from a linear isobaric equation of state of the total mass 
density p = p x + P 2 , 



p = p Q [l-a(T- 



m-a ) 



(2.5) 



for small deviations from the mean values T and £7 . For 
the room temperature solutions of ethanol in water that 
are commonly used in experiments a and (3 are positive 
[3]. 

In the continuity equation (2.1), the fluid has been 
assumed to be incompressible, i.e., the mass density p 
is constant, and the mass current is proportional to the 
divergence-free velocity field u. The reduced heat current 
Q in Oberbeck-Boussinesq approximation consists of the 
convective part uST and the diffusive part —'VST. In 
the reduced concentration current J, the diffusive part 



is —LV(SC — ipST). Here, we suppress the convective 
transport u F of the mean quantity F , since it drops 
out in the balance equations. The Lewis number L is the 
ratio of the concentration diffusivity D to the thermal 
diffusivity k, and the Prandtl number a is the ratio of 
the momentum diffusivity v and k: 



T D 

L = — ; a = — . 

K K 



(2.6) 



For room temperatures (10°C — 40°C), the Prandtl num- 
ber of ethanol- water mixtures lies between 5 and 20 [3] , 
while for normal fluid Helium it is ten or more times 
smaller. The Lewis numbers are around 0.01 for these 
mixtures. 

In experiments as well as in our simulations, the ma- 
terial parameters L and a that characterize a particular 
mixture are fixed, while the Rayleigh number R and the 
separation ratio %p, 
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(2.7) 



are considered to be control parameters that can be var- 
ied independently. R measures the externally imposed 
thermal stress. The Soret coupling ip between temper- 
ature and concentration, into which enters the thermo 
diffusion ratio kr of the mixture, reflects the influence 
of temperature gradients on the concentration field, ip 
can be positive or negative [3], depending on T and Cq. 
For room temperature ethanol-water mixtures, ^-values 
between about —0.5 and +0.2 can be easily realized ex- 
perimentally [3]. The Dufour cross-coupling term in the 
heat balance equation can be ignored for binary liquids, 
but it has to be kept for gas mixtures [84,65]. 

The buoyancy force (p — p Q )g due to density deviations 
from the mean is the driving mechanism for convective 
motion. It enters into the momentum balance (2.4) via 
the buoyancy term B which follows from (2.5) after scal- 
ing. This is the only place where a density variation is 
considered in Oberbeck-Boussinesq approximation. 

Taking the divergence of the Navier Stokes equa- 
tion (2.4), one gets, via the continuity equation, a Poisson 
equation 



vV 



-V[V(u:u)] + aRd z (ST + SC) (2.8) 



for the pressure p. The Poisson equation substitutes the 
continuity equation and builds, together with (2.2-2.4), 
a complete set of equations for the fields u, ST, SC, and 
P- 



B. Boundary conditions 

The horizontal boundaries of the layer that we shall 
call plates for shortness are at z = 0,1. The lateral 
boundaries are at x = 0, T. The plates are rigid for 
the fluid, 



u = w = for z = 0, 1 , 



(2.9) 



D. Conductive state 



and perfect heat conductors, so that the temperature of 
the fluid at z = 0, 1 is constant and laterally homoge- 
neous: 

ST = — at 2 = , and 
2 

6T = ~ at 2 = 1. (2.10) 

This is a good approximation for copper and sapphire 
plates that are used in most experiments, since their heat 
conductivity is very large compared to that of the fluid. 
The plates are impermeable to the fluid, so there is no 
concentration current through the plates: 

Je z = or 

d z SC = ipd z ST at 2 = 0,1. (2.11) 

These NSI (no-slip impermeable) boundary conditions 
have to be contrasted with the idealized FSP (free-slip 
permeable) boundary conditions that allow expansion of 
the fields in trigonometric functions and that are often 
used in theoretical approaches. However, since the latter 
allow vertical concentration transport through the plates, 
they are likely to misrepresent the delicate concentration 
balance [46] that is involved, e.g., in explaining structure, 
dynamics, and stability of LTW states. 

In lateral direction, all fields F = u, w, T, C, p are 
periodic: 

F(x,z;t) =F(x + T,z;t), (2.12) 

with a given lateral periodicity length T. Since the pres- 
sure p is determined via the Poisson equation (2.8) by 
u, ST, SC, we do not need boundary conditions for it. 

Most of the nonlinear states determined in this paper 
were obtained in a periodicity interval of T = A = 2. 
Thus the wave number is fixed at k = tt, which is close 
to those observed in experiments. 

C. Numerical method 

To integrate the partial differential equations, we used 
a modification of the SOLA code that is based on the 
MAC method [85-88]. This is a finite-difference method 
of second order in space on staggered grids for the dif- 
ferent fields, with a uniform spatial resolution Ax = 
Az = 0.05 and with an explicit first-order Euler step in 
time. The pressure field was iteratively calculated from 
the Poisson equation (2.8) using the artificial viscosity 
method [88]. Most calculations were done on a CRAY. 
For more details see [49]. 



The equations of sec. II. A together with the boundary 
conditions of sec. II. B have a stationary, laterally homo- 
geneous solution, for all parameters 

^■cond — , S T cond — — Z , 

SC cond = J, Q - z^j , (2.13) 

which is globally stable for small thermal stress. It de- 
scribes vertical heat diffusion through the layer without 
convective motion. The heat current is Q CO nd = e z • The 
constant temperature gradient induces a vertical concen- 
tration gradient (given by —ip) such that the concentra- 
tion current vanishes, 3 CO nd = , everywhere - not only 
at the plates. The concentration stratification causes in 
the buoyancy 

B cond = aR (1 + 0) Q- 2 )^ ( 2 - 14 ) 

a modification by the factor (1 + ip) relative to the ther- 
mal part. So, depending on the sign of the Soret cou- 
pling, the Soret effect enhances or depresses the ther- 
mal part of the buoyancy. The pressure is p CO nd = 

-\oR (1 + VO (\ -zf. 

To describe convection in the following, we sometimes 
use the deviations of the fields from the conductive state, 

6 = ST- ST cond ; c = SC - SC cond , (2.15) 

and similarly for the buoyancy force 

h = B-B cond = be z ; b = 6 + c. (2.16) 

E. Soret coupling 

The Soret effect causes a linear coupling between the 
degrees of freedom of the concentration field on the one 
hand and those of the temperature and also of the ve- 
locity field on the other hand. It is this coupling to the 
concentration field that causes the additional richness of 
convective dynamics and structures in binary mixtures in 
comparison to pure-fluid convection. On the one hand, 
the concentration field is mixed and adverted passively 
by the nonlinear uSC coupling. But, on the other hand, 
it enters via the Soret effect linearly into the buoyancy 
force that drives convection. Thus, the structural dynam- 
ics of the concentration field actively feeds back into the 
driving mechanism for convection. Without Soret cou- 
pling, -0 = 0, the temperature gradient of the conductive 
state does not induce a vertical concentration gradient 
- any concentration deviation from the mean diffuses 



away. And since convection does not create new con- 
centration fluctuations, the buoyancy force is not modi- 
fied. So, in the limiting case of zero coupling ip = 0, the 
equations effectively describe a pure fluid. For ip ^ , 
however, the externally-imposed temperature stress sus- 
tains — against the mixing action of convection — via 
the Soret effect concentration variations which in turn 
alter the driving force for convection. 



r osc intersect in Fig. la at the codimension-two value 
ipCTP = —3.5 • 10~ 5 with slightly different critical wave 
numbers and a small Hopf frequency. For a more de- 
tailed discussion of the codimension-two point see, e.g., 
[66-69,72,73]. 

The linear quantities of Fig. 1 do not change much as 
long as a remains significantly bigger and L significantly 
smaller than 1. For a more detailed discussion of the L, 
a dependence, cf. sec. VI. 



F. Bifurcations out of the conductive state 

The stability properties of the conductive state against 
infinitesimal convective perturbations have been dis- 
cussed in detail in the literature [66-69]. For ethanol- 
water parameters, L = 0.01 and a = 10, the results are 
summarized in Fig. 1 as a function of ip for the experi- 
mentally accessible 0-range. 

Here and in the rest of this paper we use the scaled 
Rayleigh number 



R_ 

m 



(2.17) 



that is reduced by the critical Rayleigh number R° c for 
onset of convection in a pure fluid with the critical wave 
number fc° . The analytical values are R° c = 1707.762 
and = 3.11632. They are used in Fig. 1 for the re- 
duced critical values of r and of k = k/k® . In the finite- 
differences approximation of the field equations that is 
used in our MAC algorithm, however, the threshold for 
onset of pure fluid convection [89] lies at 



R° = 1686(±0.2%) 



for 



Ax 



Az = —. 
20 



(2.18) 



So, when presenting numerical results obtained with 
the above discretization, we scale Rayleigh numbers by 
(2.18). 

The lines r s t a t and r osc in Fig. la are bifurcation 
thresholds where convective solutions branch out of the 
conductive state, and Fig. lb shows the critical reduced 
wave numbers there. At the full and dash-dotted line 
in Fig. la, a SOC solution branches off the ground 
state. For ip above (below) the tricritical value ipgoc = 
-0(1O~ 7 ) [69], the SOC bifurcation is forwards (back- 
wards). However, below ipsoc = ~ + L) [36,65] 
where r sta t diverges, the lower SOC solution branch is 
disconnected for positive r from the conductive solution 
[71]. 

The dashed line in Fig. la denotes a Hopf bifur- 
cation threshold where symmetry-degenerated left and 
right TW solutions branch out of the conductive state. 
The Hopf oscillation frequency at onset, u>h, is shown 
in Fig. lc. It varies as uj 2 h « -449 -0/(1 + ip + l/cr) 
[27,18]. For ip above (below) the tricritical value ip^w = 
-O(10~ 4 ) [69,70] the TW bifurcation is forwards (back- 
wards). The critical bifurcation thresholds r s t a t and 



G. Order parameters 

To characterize the convective solutions we use differ- 
ent order parameters, (i) The maximal vertical velocity 
w max directly measures the convective amplitude, (ii) 
The Nusselt number 



N 



dxQ z 



(2.19) 



is the total vertical heat current through the fluid layer, 

r r 

/ dxQ z , reduced by the conductive part, J dxQ con d = 

o o 

T . In our scaling, Q CO nd = 1 • We evaluate TV at z = 

and in addition at z = 1 to check that our code does 

not violate energy conservation. The reduced vertical 

heat current carried by convection alone, N — l, measures 

squared convective amplitudes of velocity or temperature 

fields. 

(iii) Since w max and N — 1 do not characterize the 
concentration field in the nonlinear convective states ap- 
propriately, we introduce a quantity that characterizes 
the magnitude of concentration variations. To that end, 
we use the mixing parameter 



M 



^(8C*)/6C, 



2 

cond 



(2.20) 



M is the variance of the concentration field reduced by its 
value in the conductive state. Here brackets and overbar 
imply lateral and vertical averages, i.e., a spatial aver- 
age over the whole fluid volume. In a perfectly mixed 
mixture where all concentration deviations SC from the 
mean vanish, M would be zero. On the other hand, in 
the conductive state with the Soret-induced concentra- 
tion gradient, M is defined to be 1. So 1 — M is an 
order parameter for the convective state that is zero in 
the ground state and approaches +1 for convection with 
strong mixing properties. 

(iv) TW states are characterized by their oscillation 
frequency to. 



III. NEGATIVE SORET COUPLING 

In this section, we discuss the structure and bifurca- 
tion behavior of extended convective states in mixtures 



with buoyancy- reducing negative Soret effect. We start 
with a typical bifurcation diagram for an ethanol-water 
mixture. We discuss the structure of the fields of the 
TW and SOC states in this diagram, including symme- 
try properties, lateral Fourier decomposition, and of time 
averaged currents of mass, heat, and concentration. The 
variation of bifurcation behavior and field structure with 
separation ratio, for the case of a buoyancy-enhancing 
positive Soret effect, and the influence of changing the 
Lewis and Prandtl number on the bifurcation structure 
are presented in Sections IV, V, and VI, respectively. 



A. Bifurcation behavior 

Here we discuss bifurcation diagrams of extended 
states in a mixture with parameters L = 0.01, a = 
10, ip = —0.25. Fig. 2 shows the r-dependence of fre- 
quency w, mixing parameter M, and Nusselt number 
N — 1 of the convective states with wavelength A = 2. 
These diagrams display in a representative way the sub- 
critical bifurcation topology that is typical for sufncently 
negative Soret coupling ip. The dotted line gives for com- 
parison the Nusselt number of the SOC states in a pure 
fluid. The conductive state loses its stability at r osc . Just 
above r osc , the system does not saturate in a state with 
small convective amplitude as in a pure fluid; rather, a 
nonlinear TW state of large amplitude is stable. This 
TW state has a frequency of only about jq^h- Reduc- 
ing the drive r quasistatically, the TW frequency con- 
tinuously increases up to about \loh at the saddle r^ w . 
Below the saddle, there is no longer a TW branch. So, 
with further reduction of the drive below r^ w , the sys- 
tem undergoes a transition to the conductive state. The 
arrows in Fig. 2c indicate the hysteretic character of the 
bifurcation. The reason for the associated discontinuous 
onset of convection is a nonlinear feedback: For ip < 0, 
the Soret-induced conductive concentration distribution 
weakens the buoyancy. Convection, on the other hand, 
redistributes the alcohol more evenly, thereby reducing 
the adverse Soret effect and thus increasing the buoyancy 
much more strongly than the flattening of the vertical 
temperature profile in the bulk decreases the buoyancy. 
The increased buoyancy in turn strengthens convection, 
which again amplifies the buoyancy. 

With increasing r, the TW frequency decreases un- 
til there is a continuous transition, at r* , to a TW of 
frequency zero, i.e., a SOC state. Approaching r* , the 
TW frequency decreases as (r* — r) = . From a fit u(r) = 
a (r* — r) h with parameters a, r* , and b, we got from our 
dataa = 0.4tt±15%, r*-l = 0.65 ±0.5%, b = 0.5 ±10%. 
Ohlsen et al. [11] obtained a slightly different exponent 
b = 0.58 from their experiments. Bensimon et al. [51] got 
b = 0.5 from their concentration-boundary-layer model. 
For larger r, we obtained stable SOC states. The transi- 
tion at r* between the upper TW branch and upper SOC 
branch (squares in Fig. 2) is smooth and non-hysteretic. 



Viewed from the SOC branch above r* , we have, with 
decreasing r, a symmetry-breaking forward bifurcation 
to left- or right-traveling waves, the symmetry of the 
latter being reduced relative to that of the former (cf. 
sec. III.C). A similar bifurcation with u> as order pa- 
rameter has been discussed in another connection, cf. eg. 
[90,91]. 

It is possible to stabilize the SOC branch below r* 
(triangles in Fig. 2) down to the SOC saddle r s soc by 
preventing this symmetry breaking. In most cases, we 
enforced the stability of the otherwise unstable SOC's by 
hindering the traveling of the rolls with the phase pinning 
boundary condition u(x, z; t) = at x = 0, T in addition 
to the periodic boundary conditions in x direction. This 
means that we do not allow lateral convective transport 
through the line x = 0. So, in our short system of T = 2, 
the roll pair cannot travel. In somewhat longer systems, 
this boundary condition is not sufficient to prevent TW's 
in the bulk, as has been seen in simulations of a system of 
length r = 10 [49] and similarly in experiments [4,28-30]. 

The Nusselt numbers of TW and SOC states typically 
lie just below the values for a pure fluid. Only near 
the saddles are there stronger deviations. The stabi- 
lized SOC's transport more heat vertically than the corre- 
sponding TW states for the same parameters. We should 
mention here that the conductive state is also with phase 
pinning boundary conditions stable only up to r osc . A 
transient standing wave of small amplitude grows above 
r osc . However, a final nonlinear standing- wave state is 
not stabilized by the phase pinning boundary conditions. 
Instead, the system ends in the SOC state. 

If we start from a stabilized SOC state (triangles in 
Fig. 2) below r* and then release the phase pinning, the 
SOC state becomes unstable. The convection rolls start 
to travel, and with growing phase velocity the Nusselt 
number decreases. Finally, for r s TW < r < r* , the system 
ends in the stable TW state. For r < r^ w , on the other 
hand, there is no stable convective state. Below r^ w , the 
convection amplitude decreases, while the phase velocity 
increases, until finally the system approaches the con- 
ductive state as a linear TW of exponentially-decreasing 
amplitude. 

That our SOC and TW states are closer to SOC states 
in a pure fluid than to the conductive state of the mixture 
is reflected not only in the Nusselt number, but also in the 
mixing parameter M (2.20) (cf. Fig. 2b) of the concentra- 
tion field. M is normalized to 1 in the conductive state 
with the Soret-induced vertical concentration variation. 
Convective mixing dramatically reduces this concentra- 
tion variation to M ss 0.1, with a slight increase at the 
saddle. There is a stronger variation of the concentration 
field in the TW's. There, M varies roughly linearly with 
frequency to, for not too small to, and reaches its maxi- 
mum (« 0.35) at the saddle r^ w . Thus, the convective 
homogeneization of the concentration, measured in terms 
of M, is more efficient by a factor of 3.5 in SOC's than in 
a TW at the saddle. A similarly-defined mixing parame- 
ter of the temperature field is reduced to only 0.9 for the 



strongest convective state in Fig. 2. 

The bifurcation diagrams in Fig. 2 were obtained in a 
system with periodicity length T = A = 2. Thus, some 
of the instability mechanisms, e.g., the Eckhaus instabil- 
ity [31-33] and more complicated dynamical behaviour 
[10,34,48,35,49,1,46, and references cited therein] occur- 
ring in large systems are suppressed here. 



B. Structure of the fields 

Here and in the next sections we discuss how the fields 
of the nonlinear convective states change along the upper 
bifurcation branch. To that end we display in Fig. 3 three 
TW states and one SOC state contained in the bifurca- 
tion diagrams Fig. 2. Their symmetry properties will be 
discussed in sec. III.C, while vertical profiles and a lateral 
Fourier analysis are presented in Fig. 4 and sec. III.D. 

Let us first discuss the fast TW close to the saddle in 
Fig. 3a. This TW travels with phase velocity v p = tu/k 
to the right. Since all its constituent fields have the form 

SF(x,z;t) = SF(x - v p t,z;0) , (3.1) 

it is a stationary structure with the given spatial period- 
icity 

SF(x + X,z;t) = SF(x,z;t) (3.2) 

with A = 2. The region -1.275 < x < 1.225 of Fig. 3a 
covers two and a half convection rolls, and the maxi- 
mum of the vertical velocity is at x = 0. The first row 
shows the colour-coded temperature field ST together 
with arrows representing the velocity field u. Close to 
the maximal upflow at i = 0, the isotherms are bent 
upwards, because warm fluid is transported upwards 
from the warmer lower plate there. Correspondingly, 
the isotherms are bent downwards close to the maximal 
downflow at x = ±1. Both velocity and temperature 
fields look like those in a pure fluid. But here in the 
mixture, the temperature field is slightly phase shifted 
relative to the vertical velocity field. The phase shift can 
be seen most clearly by comparing the lateral wave pro- 
files at z = 0.5 in the fourth row of Fig. 3. They show 
that the temperature wave ST (thin line with triangles) 
lags behind the wave of the vertical velocity w (thin line) . 
The phase shift between w and T grows monotonically 
with frequency [57], roughly like arctan(w/(fc 2 -|-7r 2 )) [55], 
and is largest for linear TWs [68] at the Hopf bifurcation 
threshold r osc . The temperature wave is not planar: in 
fact, the lines of constant phase of the first lateral Fourier 
modes of all TW fields are curved [57] to a varying degree 
instead of being straight, vertical lines. 

The third row contains the colour-coded concentration 
field S C together with the streamlines of the velocity field 
in a frame comoving with the TW. While the velocity 
and temperature fields are laterally nearly harmonic, the 
lateral concentration profile has a trapezoidal, strongly 



nonharmonic shape — see the thick line with squares for 
SC(x,z = 0.5) in the fourth row of Fig. 3a. For a TW 
propagating to the right, SC has a positive plateau with 
a concentration surplus (blue-green area) in the region of 
the right turning convection rolls and correspondingly a 
negative plateau with a concentration deficit (red-orange 
area) in the region of the left turning rolls. Between the 
two concentration plateaus there is an almost linear con- 
centration variation in horizontal direction. Also, from a 
plateau to a plate there is a steep, almost linear variation 
of SC (cf. Fig. 4c). 

To explain this structure of the TWs concentration 
field self-consistently, we view the fluid layer from the 
frame of reference £ that comoves with the phase veloc- 
ity v p e x of the TW. In this frame, all fields SF of the 
stationary TW 

u(x, z) = u(x — v p t, z; t) — v p e x (3-3) 
Sf(x,z) = ST(x-v p t,z;t) (3.4) 
SC(x,z) = SC(x — v p t,z;t) (3-5) 

are time- independent, with 

x = x — v p t . (3-6) 
From the conservation law 

St SC + V • J = (3.7) 
and d t SC = 0, it follows that the concentration current 
3 = uSC-LV (SC-*I)ST} (3.8) 

is divergence-free. Then the continuity equation (2.1) 
implies 

(u- V) SC = L V 2 (SC - i) ST ) . (3.9) 

Since the Lewis number L is small, we can set the dif- 
fusive term on the r.h.s. of (3.9) to zero for the present 
discussion. Thus, there is no gradient of SC parallel to u; 
i.e., the isolines of SC lie practically parallel to u. Fur- 
thermore, if we define the velocity streamfunction $ in 
£ via 

-Vx$ e3/ = (d z $,0,-d x $) = (u,0,w), (3.10) 

then the isolines of <3>, i.e., the streamlines of u along 
which the fluid moves in S, are almost identical with the 
isolines of the concentration. That this holds quite well 
can be seen in Fig. 3 by comparing the colour-coded con- 
centration distributions with the dashed superimposed 
streamlines of u. So, in order to understand the struc- 
ture of the concentration field, we first investigate the 
form of the streamlines. The latter can be most easily 
understood within a single-mode representation [55], 



$ = — — w(z) sin(fca;) — v p z . (3-11) 

fx, 

This streamfunction captures the characteristic proper- 
ties of the velocity field being basically harmonic in the 
lateral direction with a vertical profile w(z). The second 
contribution, — v p z, to (3.11), that arises in £ from the 
phase velocity of the TW, has two effects: (i) it shifts 
the closed isolines of the ordinary first contribution to 
(3.11) — corresponding to stationary roll convection — 
alternatingly towards the top and bottom plates, and 
(ii) it generates open streamlines. These meander be- 
tween and around the roll-like regions of closed stream- 
lines (Fig. 3a), whereas in an SOC (v p = 0), the only 
open streamlines are the vertical separation lines be- 
tween the oppositely turning rolls joining top and bot- 
tom plate. Consequently, in a right-propagating TW, (i), 
the blue-green, concentration rich (red-orange, concen- 
tration poor) regions of closed streamlines for the right 
(left) turning fluid domains in £ are displaced towards 
the upper cold (lower warm) plate, where the Soret ef- 
fect has caused a concentration surplus (deficiency) for 
ip < 0. In addition, (ii), open streamlines separate the 
blue-green (red-orange) regions from the bottom low- 
concentration (top high-concentration) boundary layer. 
Therefore, the blue-green (red-orange) boundary layers 
feed high concentration (low concentration) only into the 
respective rolls. See the downwards blue (upwards red) 
jet-like structures of the concentration field pointing into 
the two regions of closed streamlines in Fig. 3a. Within 
the regions of closed streamlines, the concentration field 
behaves like a passive scalar, mixing up diffusively [78] 
as the fluid is advected around in £ along the closed 
streamlines. Since a volume element in the area of closed 
streamlines moves on a closed orbit in S, the fluid in 
such an area is caught within it if we ignore diffusion. 
Thus, in the laboratory frame, this entire fluid volume 
is transported with the phase velocity v p e x of the TW 
to the right. So the concentration wave moves about j 
wavelength ahead of the w-wave (see the profiles of the 
fields in the 4th and 5th row of Fig. 3). 

The above explanation of the concentration field struc- 
ture is based on the smallness of L and on the fact that 
SC shows some of the characteristic behaviour of a pas- 
sively convected field. If one tried to explain the structure 
of the temperature field along similar lines, starting from 
the divergence freeness, V • Q = 0, of the heat current 

Q = uST -VST (3.12) 

in S, one would see that in this case (for |u| of the order 
of 10) the diffusive term in 

(u- V) ST - V 2 6T = (3.13) 

cannot be assumed to be small compared to the convec- 
tive term as in the concentration case. 



Let us now discuss the TW variation with decreas- 
ing frequency, i.e., increasing r. Since the phase 
shift between velocity and temperature wave varies like 
arctan(w/(fc 2 + 7r 2 )) [55,57], it vanishes almost linearly 
with uatr*. From the streamfunction $ (3.11), one fur- 
thermore infers that, with decreasing phase velocity v p 
of the TW, (i) the alternating up- and downwards shift 
of the closed streamline regions and (ii) the transverse 
extension of open streamlines of u between the rolls as 
well as between a roll and a plate decreases. Both ef- 
fects lead to a reduction of the concentration contrast 
between the differently displaced and rotating regions 
of closed streamlines: decreasing v p makes the stream- 
line structures more and more mirror symmetric around 
x = 0. Secondly, the feeding of concentration surplus 
(deficiency) from the top (bottom) boundary layer into 
the different roll regions becomes more symmetric. 

Let us now try to understand the structure of the slow 
TW's and the SOC states. The picture of mixing of the 
TW concentration field within a region of closed stream- 
lines of u was based on the assumption that the convec- 
tive term in the concentration balance dominates. For 
small phase velocities, and especially for SOC states, this 
picture has to be extended, since these closed streamlines 
cover an entire convection roll (see third row in Fig. 3d). 
Since the velocity field vanishes at the plates, concen- 
tration diffusion has "long enough time" to work before 
a fluid volume is transported away from the plates by 
the velocity field, and a Soret-induced concentration gra- 
dient develops at the plates. The turnover time of a 
fluid element that starts very close to a plate on a closed 
streamline is indeed much bigger than the vertical con- 
centration diffusion time j- . Since the spatial area where 
the turnover time is at least comparable to is small, 
the concentration boundary layer at the plates is small. 
Since streamlines through the top layer also pass through 
the bottom boundary layer, a fluid element on such a 
streamline experiences a strong change of its concentra- 
tion content. Let us now follow a volume element starting 
close to the upper plate, say, at x = 0.5 with positive SC 
and negative ST. First, it moves slowly with decreas- 
ing velocity to the right roll boundary while retaining 
its temperature and concentration. After that, it moves 
slowly away from the upper plate, since there w is just 
proportional to the square of the distance from the plate. 
Then it is transported very quickly downwards through 
the middle of the layer, where w max is of the order of 10. 
Since heat diffusion works quickly, the element warms up 
rapidly, and the isotherms there are bent only moder- 
ately downwards, as in a pure fluid. However, with con- 
centration diffusion working very slowly on the timescale 
of -|- = 100 thermal diffusion times, the fluid element re- 
tains positive concentration down to the vicinity of the 
lower plate. This causes the little peak in the lateral pro- 
file of SC in Fig. 3d at the position x = ±l,z = 0.5 of 
maximal downflow. In the vicinity of the lower plate, the 
volume element is transported very slowly further to the 



lower plate and then to the left. In this part of the cycle 
concentration diffusion again has enough time to work. 
Surplus concentration is now transported diffusively - 
mostly upwards - out of the volume element, and the 
concentration content of the fluid volume becomes nega- 
tive. All these effects together generate the typical plume 
structure of the SOC concentration field with the shafts 
of the plumes centered at the positions of maximal up- 
and downflow. Such structures have also been seen in 
simulations of thermohaline convection [74,75]. The con- 
centration plumes have rather narrow vertical extension 
and modify the concentration field in the interior of the 
rolls only very slightly. Similar plumes occur in the tem- 
perature field in very high Rayleigh number convection in 
a pure fluid, as a consequence of the much faster convec- 
tive transport of heat compared to heat diffusion, where 
Veronis [92] has also seen them. 

While the plumes, or concentration jets, in the SOC 
are mirror symmetric around x = 0, this is not the case 
in TW's. In particular, in the fast TW of Fig. 3a, the 
concentration jets alternatingly feed only the left or right 
turning rolls. There, the plumes are strongly bent into 
the respective rolls, and the resulting fine structure of 
concentration enhancement (depletion) causes the peak 
(dip) at x ~ 1 (a; ~ 0). Since the turnover time also di- 
verges [55,79] on the separatrix of the closed streamlines 
in a TW, this peak (dip) fine structure on top of the 
smooth background of the interior of the rolls thus has 
the same origin in TW's as in SOC's. 



C. Symmetry properties 

The concentration field most clearly shows that SOC 
and TW states discussed here are strongly nonlinear 
structures. Nevertheless, they have, in addition to their 
stationarity (3.1) and spatial periodicity (3.2), an addi- 
tional mirror-glide symmetry 

SF(x,z;t) = ±SF^x + ^,l-z;tj (3.14) 

with + for u, p and — for w, ST, SC and B. This is a 
combination of a translation by half a wavelength in the 
x direction and a reflection through the z = \ plane. 
Besides (3.1) and (3.2), this is the most elementary, i.e., 
longest persistent symmetry of SOC's and TW's. As far 
as we know, the mirror-glide symmetry (3.14) has first 
been discussed in ref. [54]. Veronis found that Fourier 
modes breaking this symmetry decayed to zero for ther- 
mohaline convection with fixed spatial phase (i.e. no TW 
was possible) [74] and for SOC states in one-component 
fluids [92]. But the symmetry was not explicitly stated. 
The mirror-glide symmetry of SOC rolls was recently ver- 
ified numerically in convection of one-component fluids 
[89] as well. 

The symmetry (3.14) can easily be seen in the light 
intensity distribution 



I(x,z)=AV 2 [SC(x,z) +bST(x,z)] . (3.15) 

that one would observe in side-view shadowgraphs ob- 
tained by shining light horizontally through the fluid 
layer along the roll axes [93,94,3,14,15]. A is a con- 
stant, and b = —0.919 with our parameters for a real 
mixture [14,3]. Fig. 3 shows this light intensity distribu- 
tion in grey scale. Note that its characteristic structures 
are dominated by the large Laplacians of the concentra- 
tion field. This incidentally shows that the shadowgraph 
method also gives information about the concentration 
field [14,15]. 

Applying the operation (3.14) twice yields the pure 
spatial translation symmetry (3.2) by one wavelength. 
On the other hand, if one combines (3.14) with the sta- 
tionarity (3.1), one gets a space-time symmetry 

SF{x,z;t) = ±Sf(x, l-z;t+^j , (3.16) 

with r = 27r/w for the TW states. Here the transla- 
tion in (3.14) by half a wavelength is replaced by a time 
translation by half a period. This space-time symmetry 
was mentioned in another connection by Weiss [95] . Note 
that (3.16) does not imply a horizontal mirror symmetry 
at the midplane for SOC states - which they do not have 
anyhow - since r diverges for u = 0. But in a system 
S, that moves with constant velocity V e x relative to the 
laboratory system, the spatial symmetry (3.1) becomes 
a space time symmetry 

SF(x,z;t) = ±SF (x, 1 -z;t + ^ (3.17) 

with f = X/V. This also holds for TW's with V being 
the velocity of £ relative to the rest system £ of the TW. 

The only additional symmetry of SOC's compared to 
TW's is the mirror symmetry 

8F{x, z) = ±SF(-x, z) (3.18) 

with — for u and + for w, ST, SC, p, and B, where x = 
is the position of a roll boundary. All other symmetries of 
the SOC states follow from (3.1), (3.14), (3.18). For ex- 
ample, the point mirror symmetry around the roll center 
r = — 7 = — 

SF(x, z) = ±SF Q - x, 1 - zj , (3.19) 

with + for p and — for u, w, ST, SC, and B, is a com- 
bination of (3.14) and (3.18). Finally, we mention that 
linear SOC and TW transients have no further symmetry 
besides those discussed here. 



D. Vertical profiles and lateral Fourier 
decomposition 

In this section, we continue to discuss how the struc- 
ture of convection changes with frequency and Rayleigh 



number. To that end, vertical profiles of the concen- 
tration field at different lateral positions are shown in 
Fig. 4c, and the profiles of lateral Fourier modes of the 
temperature and concentration field are shown in Fig. 4a 
and Fig. 4b, respectively, for the three TW states and 
the SOC state in Fig. 3. The Fourier decomposition of 
the fields for these states is 

oo 

SF(x,z;t) = SF (z) + ReY^ SF n (z;t) e ink x (3.20) 

n=l 

with k = Because of (3.1), the moduli of the Fourier 
coefficients, 

6F n (z;t) = \6F n (z)\e- i( ""W, (3.21) 

are time independent, while in a TW the phases <pp n vary 
linearly in time. Fig. 4a shows the vertical variation of 
the Fourier modes |0 n (2)| = |<^ra(z)| of the temperature 
field for n = 1, 2, 3 and of O (z) = Sf (z) - {\-z), 
i.e., the laterally averaged deviation from the conductive 
profile. 

Due to nonlinear effects, the profile of \6\\ flattens as 
convection intensifies with growing r. The nonlinear ly 
stimulated higher harmonics \6^\, \8s \ grow with r but are 
weak over the r-region considered here. So, as in a pure 
fluid [83,89], the temperature field of SOC's and TW's 
can be accurately approximated by two lateral Fourier 
modes 

ST(x,z;t) ~ ST (z) + |iTi(z)| cos (kx - ip Tl (z; t)) . 

(3.22) 

For a TW, the phase profile 

<p Tl (*;t) =<p Tl (z;t = 0) + kv p t (3.23) 

is ^-dependent [54], while, in a SOC, tp^ is constant in 
z and t. The vertical derivative of the mean deviation 9q 
from the conductive temperature field at z = 0, 1 gives 
the convective heat current. All in all, the temperature 
field is similar to that in a pure fluid [89]. This also holds 
for velocity modes which are not displayed in Fig. 4. 

The behavior of the concentration field is drastically 
different. It shows four characteristic features: (i) con- 
vective reduction of the Soret-induced difference between 
top and bottom plates, (ii) homogeneous distribution 
in the roll-like regions of closed streamlines of Fig. 3, 
(iii) linear boundary-layer variation, and (iv) plumes. Al- 
ready close to the TW saddle, the concentration differ- 
ence between the plates, C(z = 1) — C(z = 0), is convec- 
tively reduced to about | of the value for the conductive 
profile (dotted line in Fig. 4c). This ratio is further re- 
duced to | for the SOC state in the right column. The 
reason is that impermeability [27] imposes a condition 
only on the vertical derivative of the concentration field 
at the plates, while the temperature is fixed via a Dirich- 
let condition at the plates. A similar reduction of the 



concentration difference between the plates has also been 
found in the Galerkin model of Linz et al. [71,55] for im- 
permeable boundary conditions. The concentration cur- 
rent at the plates, 

3(x,z;t) = -Ld x SC(x,z;t)e x for z = 0, 1 , (3.24) 

is horizontal and small, but in general nonzero. 

The homogeneous distribution in the roll-like regions 
in Fig. 3 shows up in Fig. 4c as plateaus of the verti- 
cal profiles of SC. Note that the profiles A, B are re- 
lated to profiles A', B' via the symmetry (3.14). Curve 
A (A') refers to a section through a right (left) turning 
roll at xa = 0.35 (xa> = —0.65) close to the roll center 
at \ (-i). Profile B at x B = 0.85 (B' at x B > = -0.15) is 
slightly further away from the center, so that the vertical 
extension of the plateau region is slightly smaller there. 
Since xb {xb 1 ) is close to a downflow (upflow) jet from 
the concentration rich top (poor bottom) plate, curve B 
(B') lies in the upper (lower) part of the layer above (be- 
low) curve A (A'). The high concentration blue (low con- 
centration red) plume that is best visible in the TW of 
Fig. 3a leads to a concentration enhancement (depletion) 
in the low- z (high- z) part of the vertical profile A (A'). 
This consequence of the previously-described plumes can 
also be seen in the profiles of the slow TW's and SOC of 
Fig. 4c. 

Since the convective mixing of concentration is very 
effective, the laterally-averaged concentration (6C)(z) is 
nearly zero over a wide range in the bulk of the layer. 
For our stationary TW's and SOC's, (•) denotes a lateral 
average over a wavelength 

(f)(z) = jf*f(x,z;t)dx, (3.25) 

which for TW's is equivalent to a time average over a 
period 

</>(*) = - I' J(x,z;t)dt . (3.26) 

T JO 

The impermeable boundaries cause a vertical concentra- 
tion gradient 

d z (SC) (z) = -tpN for z = 0, 1 (3.27) 

at the plates. Although it is larger than the gradient, 
d z C C ond = —tp, m the conductive state (dotted line in 
Fig. 4c), the overall concentration difference between the 
plates is strongly reduced. In this way, the density differ- 
ence between the plates and with it the buoyancy force in 
the fluid has increased, as compared with the conductive 
state. For Rayleigh numbers larger than those presented 
here, mixing effects can also be seen in the temperature 
field. Together with those of the concentration, they lead 
to a mixing in the mass density, so that there is no ver- 
tical gradient of the mass density in part of the area of 
closed streamlines (cf. Fig. 3b of [54]). 



As a consequence of the reduction in the concentration 
difference between adjacent rolls with decreasing TW fre- 
quency, the strength of the concentration Fourier modes 
decreases, in particular in the middle of the cell. Since 
the area of the open streamlines of u also shrinks, the 
plateaus of \SC\\ broaden (thick lines in Fig. 4b). In the 
SOC state, SC nearly vanishes in the middle of the cell, as 
can already be inferred from the smallness of the mixing 
parameter M (Fig. 2b). Only the boundary layers be- 
tween the rolls and the plates remain. However, the lat- 
eral variation of the concentration right at the plates and 
in its immediate vicinity — see the little orange (blue) 
dip at the origin of the upwards jet at x = (downwards 
jet at x = ±1 ) in Fig. 3d — implies that SC\ is finite 
at z = 0, 1 (cf. Fig. 4b). In fact the SOC concentration 
right at the plates obeys a Laplace equation 

(d 2 x + dl) SC(x, z) = for z = 0, 1. (3.28) 

The boundary-layer behavior of the concentration also 
shows up in the nearly linear variation of the laterally- 
averaged concentration fields (SC) (z) = SCq(z) (dashed 
line in Fig. 4b) at the plates with a gradient —ipN. 
The consequence of the plume structure caused by bend- 
ing the upwards (downwards) concentration jet laterally 
along the plate can be seen in the SOC state of Fig. 4b 
by the weak minimum (maximum) of SCq below (above) 
the boundary layer at the upper (lower) plate. Thus, 
SCo remains slightly negative (positive) until z = \. In 
these regions, the concentration gradient is inverted, giv- 
ing rise to a weak destabilizing effect in the buoyancy 
force. Inverted concentration gradients have also been 
seen in thermohaline calculations [74]. Such an inver- 
sion is again more common as an effect in the tempera- 
ture field in high Rayleigh number convection in a pure 
fluid, e.g., [92]. Note also that slow TW's (third column 
in Fig. 4b) show a weak inversion of the concentration 
gradient. For faster TW's (first and second column in 
Fig. 4b) , the large boundary-layer behaviour in the open 
streamlines of u dominates. 

The Rayleigh number variation of the Fourier modes 
at selected vertical positions is displayed in Fig. 5. For a 
pure fluid, the velocity modes |w2|, and \w%\ grow 

with r [89]. In our binary mixtures, this also holds suf- 
ficiently far away from the TW saddle. Approaching the 
latter from above, in particular |w?3| shows a marked in- 
crease (full triangles in Fig. 5a), which is combined with 
a strong, but smooth phase shift for the TW's, while 
for the enforced SOC's, wz goes through zero, combined 
with a phase jump by tt. There is no significant change of 
the behavior at the transition from the TW to the SOC 
branch at r* . In a pure fluid, the temperature modes 
\6i\, \02 1, and |#3| also grow with r. Note that the tem- 
perature is scaled with AT, so that \9i \ decreases with r 
above a value r ~ r* . At the saddle, there is only a small 
difference for the states in the mixture compared to pure 
fluid convection. However, one can see a strong increase 



in the concentration modes close to the saddle. The rea- 
son for this increase lies in the smaller velocity field at the 
saddle and with that a smaller mixing effect in the con- 
centration field. The structure of the SOC states with or 
without phase-fixing boundaries is the same on the entire 
upper branch with a small increase of the concentration 
modes close to the saddle r s soc . 

The variation of concentration Fourier modes with TW 
frequency is shown in Fig. 6. In the bulk of the layer, 
\SC\\ grows roughly linearly with u (Fig. 6a), which was 
also observed in experiments [15]. The deviation of \SC\\ 
at z = 0.25 from the bulk value at z = 0.5 for large u> 
reflects the increase of open streamlines of the velocity 
field when approaching the TW saddle. At the plates 
(see z = 0) we always have boundary-layer behavior, so 
that the concentration variation with to is small. The 
concentration wave is phase shifted somewhat more than 
j with respect to the velocity wave in the bulk of the 
layer (Fig. 6b). This is also seen in the experiments of 
Winkler and Kolodner [15]. For very small frequencies, 
there is an increase of the phase shift in the bulk to ^ for 
SOC states. At the plates (see z = 0), however, the con- 
centration wave is shifted significantly less than j to the 
velocity wave for TW states and shows no phase shift for 
SOC states. The laterally-averaged concentration field is 
nearly zero in the bulk (see Fig. 6c for z = 0.25). Thus, 
the shrinking of the area of closed streamlines and the 
corresponding smaller mixing area of the concentration 
has only a weak influence there. This smaller mixing 
effect can be seen better in the concentration contrast 
between the plates, which increases linearly with u> for 
not too small to. 



E. Mean flow and time averaged currents of heat 
and concentration 

In this section, we discuss mean transport properties 
of TW's. In a TW, mass (fluid) is transported along the 
streamlines of u in the comoving frame S. Thus, in the 
laboratory system, part of the fluid moves to the right, 
while another part moves to the left. Here, however, we 
do not consider this Lagrangian motion [8,9,57,56,55,79] 
but the time-averaged mass transport through an area S, 

(^J ds- pj = p Q J ds- (u) , (3.29) 

taking the mass density to be constant. So we discuss the 
time-averaged velocity field (u)(x, z). In a SOC state, 
this time average once again yields the velocity field. For 
a TW, temporal and spatial lateral averages are equiva- 
lent, because of (3.1), so that (u) and all other averaged 
quantities (SF) depend only on z and no longer on x. 
Because of the incompressibility of the fluid, (w) = 0. 
Thus, streamlines of (u) (Fig. 7a) run horizontally. The 
mean flow [55] 



U(z) = (u(x,z;t)) 



(3.30) 



determining the time averaged lateral mass current 
PoU(z) [54,57] is opposite to the phase velocity of the 
TW for all z. The vertical average U = U(z) dz of 
this mean flow velocity is very small compared to v p and 
to the maximal velocity in the fluid, e.g., U ~ —0.008 for 
a TW at r = 1.246, -tp = -0.25. As shown by Linz et 
al. [55], the mean flow U is related to the curvature of 
the phase, p Wl (z), of the dominant velocity mode. See 
Ref. [54] for the profile of ip Wl (z) and of U{z) [96]. A 
small curvature of p Wl is also present in linear TW's 
[66,54,68]. But there the curvature is such that the mean 
flow is parallel to the phase velocity [68]. This is compat- 
ible with the r-variation of \U(z = 0.5) | shown in Fig. 8a. 
This extremal value of the mean flow does not only go to 
zero at r* but also seems to approach zero in the neigh- 
borhood of the saddle r^ w . So we expect from the linear 
theory [54,68] that the mean flow changes sign near r^ w 
and has weakly positive values on part of the lower TW 
branch and at r osc . 

The time-averaged heat current, 



(Q) (z) = Ne z + (STu) (z)e x , 



(3.31) 



is dominated by its vertical contribution TV e z . Also, the 
mean convective heat current 

(Q - Qcond) (z) = (N — 1) e z + (ST u) (z) e x (3.32) 

is dominated by the vertical contribution. But the 
streamlines of (Q — Qcond) are bowed in the bulk of the 
fluid into the direction of the phase velocity of the TW 
(Fig. 7b): in addition to the z- independent vertical part, 
the TW generates a z- dependent mean lateral heat cur- 
rent (STu) [72,55,54,57], which is parallel to the phase 
velocity in the lower half of the fluid layer and opposite 
to the phase velocity in the upper half. Using a lateral 
Fourier decomposition of the fields, the averaged lateral 
heat current is [72,55,54] 



(Qx) (z) = (STu) 

1 °° 1 - 
= (ST) U+-lmJ2 ~ST n d z w* n 



(3.33) 



n=l 



Since U(z) and also the higher lateral harmonics of w 
and ST are small, the main contribution to (STu) comes 
from the first lateral harmonic [55,54,57,48], i.e., 



(STu) (*)-^ 



ST i 



(d z |wi|) sin(<p Tl -tp Wl ) ■ (3.34) 



So the lateral heat current is driven by the phase shift 
between the temperature and the velocity wave [72,55], 
which by the way also exists in linear TW's [54,68]. 
Since for these small frequencies we have ip^ ~ <Pw! ~ w 
[72,55,57], it is evident that the lateral heat current 
(Q x ) is nearly proportional to the TW frequency. This 



can be seen directly by comparing the maxima of (Q x ) 
shown in Fig. 8b as triangles with the TW frequency 
in Fig. 8c. The curvature of p^ (z) is not very impor- 
tant for the heat current, since it is much smaller than 
(i^Ti — <Pw 1 ) ( z = 0.5). As ip Wl (z) it is bowed on the 
upper TW branch opposite to the propagation direction 
of the TW, but to a stronger degree [54]. However, for 
linear TW's close to the bifurcation threshold the curva- 
ture of ipTi (z) like ip Wl (z) is bowed into the propagation 
direction of the TW's [68]. 

Similar as for the heat current, the phase shift between 
the concentration and velocity waves drives an averaged 
lateral concentration current [72,55,54] 



<Jx) (*) 
-1 

~ 2k 



(SCu) 

(d z |wi|) sin (<p Cl -<p Wl ) ■ 



(3.35) 



Like the lateral heat current, it is roughly proportional 
to the frequency — compare the dots in Fig. 8b which 
represent the maxima of (J x ) with the TW frequency in 



and 



Fig. 8c. But here the u> dependence arises from 

not from the phase difference: the concentration wave 
is in the bulk of the fluid about a quarter wavelength 
ahead of the w-wave for all TW's considered here (cf. 
sec. III.D). Since concentration is conserved and because 
of the impermeability of the plates, there is no averaged 
vertical concentration current, i.e., 



(J) (z) = (SCu) (z) B X , 



(3.36) 



and the streamlines of (J) (Fig. 7c) run horizontally. 
Note that this concentration current has nothing to do 
with the mean flow U. The former is in the upper half of 
the fluid layer parallel and in the lower antiparallel to the 
phase velocity [54,57]. So the lateral concentration cur- 
rent runs opposite to the lateral heat current. The vari- 
ation of ipd (z) is stronger and more complicated than 
for the temperature and velocity field since the concen- 
tration field shows strong boundary layer behavior at the 
plates, cf. Fig. 6b and [54, thick line in Fig. 2c]. For a 
more detailed comparison of linear and nonlinear TW's 
cf. [68]. 

The symmetry of the time-averaged quantities 



(SF) (z) = ± (SF) (1 - z) 



(3.37) 



follows directly from (3.14) with + for (u),(Q z ) and — 
for (Q x ), (J x )- Vertically integrating the mean currents, 
one finds that the mass current, poll e x , is very weak, the 
total heat current, 

(Q) = 7Ve z , (3.38) 

is strong, and the total concentration current vanishes: 

(J) = . (3.39) 

Some questions concerning the significance of these 
currents in extended convective states are not answered. 



The mean flow U is very small - does it play a role? 
How does the associated mass current look in a finite 
cell? Does it have an influence there? Does one have 
to consider spatial density variations [55,57]? This goes 
beyond the Oberbeck-Boussinesq approximation, which 
is the basis of our calculations. If one takes these den- 
sity variations into account in the averaged lateral mass 
current, one should do this also in the vertical mass cur- 
rent without violating mass conservation [57]. The heat 
current is a local current. But what role does the lateral 
large-scale concentration current play for TW's? Because 
of concentration conservation, it has to turn around at 
the lateral cell boundaries. Is the deformation of the 
streamlines at the boundaries more important? In any 
case, the concentration current plays an important role 
in the dynamics of LTW's [46]. 



IV. CHANGING THE SORET COUPLING 

In sec. Ill, we discussed a typical bifurcation diagram 
and the structure of convective states for ip = —0.25. 
Now we discuss how the bifurcation behavior changes 
in ethanol- water mixtures with different negative ip. We 
also show the TW and SOC structure for a small negative 
ip = -0.01. 

For experimentally accessible [3] negative separation 
ratios —0.6 < ip < —0.005, we display the bifurcation 
behavior in the tp-r plane in Fig. 9b. To resolve different 
orders of magnitude, logarithmic scales are used. For 
orientation, the bifurcation threshold of TW's, r osc , is 
included. In the ip region considered here, the conductive 
state is stable for r < r osc . Consider first the behavior 
for our reference value ip = —0.25. The r region with 
stable TW's is limited from below by the saddle ri rw 
(full dots in Fig. 9b) and from above by the transition 
to stable SOC states at r* (full squares in Fig. 9b). The 
upper SOC branch is stable only above r*. But it can be 
stabilized down to the SOC saddle r s soc (open squares 
in Fig. 9b) by using phase pinning boundary conditions 
that prevent TW's. The band limits r* ,r s TW of stable 
TW's and r osc grow as ip becomes more negative. We 
found 



All quantities r Q 



and r?p W decrease with grow- 



0.5 < 



4wW - 1 

r OS c(ip) ~ 1 



< 0.7, 



(4.1) 



i.e., the TW saddle varies with ip similarly to r osc (ip). 
On the other hand, the upper boundary r* increases 
strongly at negative ip: r*(ip = —0.25) ~ 1.65, while for 
ip = -0.295 TW's exist up to r ~ 5, and for ip = -0.3 the 
Soret coupling is so strong that TW's exist there at least 
up to r = 10 « 7r osc . For reasons of numerical accuracy 
and stability, we did not investigate larger r. This very 
strong increase of r* suggests the existence of a thresh- 
old ip below which no SOC state is stable. Except for 
this extreme increase of the region of stable TW's, there 
is no qualitative difference in the bifurcation structure 
compared to ip = —0.25. 



ing ip, and r* and r osc intersect at ip* ~ —0.09. Thus, 
upon crossing the Hopf bifurcation threshold r osc , the ini- 
tial transient growth of the fast propagating linear TW 
saturates into a nonlinear slow TW (SOC) state on the 
upper bifurcation branch for ip < ip* (ip > ip*). Fig. 9a 
shows the frequency of these nonlinear states at r osc . The 
TW frequency at r osc could be fitted with the ansatz 
^{r osc {iP)) = a \n(iP/iP*) using a = 0.25 ± 3% , iP* = 
0.091 ± 0.5%. This corrects the wrong formulation in 
footnote 1 on page 92 of [49]. Tab. I lists besides r osc 
and loh properties of the nonlinear state at r osc . With 
further increase of ip, r* approaches the saddle r s TW , so 
that for ip = —0.01 TW's exist only in a very small in- 
terval, 10~ 5 < Ar < 10~ 4 , close to the saddle. Further- 
more, the SOC saddle r s soc is only just below the TW 
saddle. But for ip = —0.005 we did not find any stable 
TW. We did not investigate smaller negative ip. Thus, 
we conclude that the end point of the TW branch r* lies 
on the upper SOC branch for ip < —0.01. For negative 
ip > —0.005, the endpoint has moved around the SOC 
saddle and lies on the lower SOC branch. So we cannot 
stabilize the SOC's further with phase pinning boundary 
conditions, and we identify the SOC saddle with the last 
stable SOC. It should be mentioned that decaying SOC's 
to the left of the stable SOC saddle eventually end as 
decaying to the conductive state as standing waves. 

The structure of TW and SOC states for small Soret 
coupling, ip = —0.01 (Fig. 10), is in principle the same as 
for the strong one, ip = —0.25 (Fig. 3). The main differ- 
ence is that the concentration variation is much smaller, 
and that there exist stable TW's and SOC's at a much 
smaller drive r. So the strength of convection and of the 
variation in the temperature field (top row of Fig. 10) is 
much smaller than for the states with the higher drive 
(top row of Fig. 3). Similarly, the convective part of the 
concentration current is not as dominant as for stronger 
drive. So the concentration boundary layer at the plates 
and between the rolls is broader (third row of Fig. 10), 
and the little dips and peaks on the concentration pro- 
files of Fig. 3 have grown here (fourth row of Fig. 10) 
into broad bumps, and the plume- like structure is much 
"softer" here than for the case ip = —0.25 with strong 
convection. Since the TW frequency is very low, the 
area of open streamlines of u is small and with it also the 
concentration difference between the rolls. So the concen- 
tration field is governed more by the plume- like structure 
than by the concentration plateaus. Since the concentra- 
tion field is weaker and smoother than for ip = —0.25, 
it contributes less to the sideview shadowgraph intensity 
(second row of Fig. 10), which hence reflects predomi- 
nantly the smooth harmonic temperature field. 



V. POSITIVE SORET COUPLING 

As mentioned in the introduction, squares and oscil- 
lations between perpendicular rolls can be observed for 
ip > 0. We cannot describe these 3d effects with our 
2d algorithm. Our calculations should be compared with 
measurements in narrow convection channels that favour 
parallel rolls. Here we show among other things that 
the structure of the concentration field of SOC states, 
with not too small amplitudes, can be explained by the 
same mechanisms for positive ip that we have discussed 
for negative ip. 

For positive ip, there is only a stationary instability 
(Fig. 1). The bifurcation threshold r s t a t becomes small 
for not too small ip, while k c drops to zero, so that the 
wavelength with largest growth rate in a real system is 
determined by its length. Here we consider ip = 0.1, 
where k c = and r stat (k c ) = 0.04216. The stabil- 
ity curve r sta t(k) is very flat, so that the bifurcation 
of SOC rolls of wavelength A = 2 occurs already at 
r s tat(k = 7r) = 0.06012 [68]. We have calculated four 
convective states in a system of length T = 20 with later- 
ally periodic boundary conditions. We have chosen this 
large T to allow for wavelength selection. As initial con- 
ditions, we took in each case a localized traveling wave 
at ip = —0.25 [46]. In each case, a SOC state of wave- 
length A = 2 finally developed (squares in Fig. 11), thus 
suggesting the existence of a strong selection mechanism. 
Only the state with the biggest r-value was calculated in 
a system of length T = 2. The first two states in Fig. 11 
are in the so called Soret region, r < 1, and the others 
in the Rayleigh regime, r > 1. Note that even the small- 
est r-value (r=0.506) lies distinctly above the bifurcation 
threshold. As expected [26], the convective heat current 
N — 1 (Fig. lib) in the Soret region is very small, while 
in the Rayleigh region it is comparable to the current in 
a pure fluid. The mixing parameter M (Fig. 11a) indi- 
cates that the fluid in the Rayleigh regime is significantly 
better mixed than in the Soret regime. But M being well 
below 1 clearly shows that the states in the Soret regime 
also have to be considered as " strongly nonlinear states" . 

Now we discuss structural variations of SOC states 
with r. Consider first the state in Fig. 12d at r = 2.026 
with the largest amplitude. Its concentration is well 
mixed in the rolls as in the SOC states at negative ip 
(cf. Fig. 3d). Also, the concentration boundary layer is 
present, as at negative ip. But here SC is positive at the 
lower plate and negative at the upper plate, since ip > 0. 
Also, the little peaks in the lateral concentration profile 
(fourth row of Fig. 12d) now lie at the position of maxi- 
mal Mpflow and symmetrically the valleys at the positions 
of maximal downflow. 

With decreasing r (see Fig. 12c, r = 1.114), the ve- 
locity field decreases, the concentration boundary layer 
grows, and the little peaks become higher and broader. 
But mixing in the rolls remains dominant down to the 
beginning of the Soret region (see Fig. 12b, r = 0.912). 



However, here the plume structure is very soft, as in the 
SOC state at ip = —0.01 in Fig. 10b, since the velocity 
field is very weak. 

With further decrease of r and thereby of the velocity 
field, concentration diffusion becomes comparable to ad- 
vection in a larger spatial area, the concentration plumes 
become very broad (Fig. 12a, r = 0.506), and the lat- 
eral concentration profile (third row of Fig. 12a) becomes 
smooth. Since convective mixing is weak, the ampli- 
tude of the lateral concentration profile at z = \ has 
increased, and the vertical concentration difference be- 
tween the plates has also increased to almost half the size 
in the quiescent conductive state. With further decrease 
of r, the nonlinear advective term in the concentration 
balance loses more and more dominance and then rele- 
vance, and the mixing parameter M increases towards 1 
at the bifurcation threshold. 

The second and last rows of Fig. 12 show the shad- 
owgraph intensity I(x,z) (3.15) with b = —0.919. The 
temperature field dominates in the Rayleigh regime, and 
the concentration field adds sharp features to it, as for 
not too weakly negative ip (Fig. 3d and ref. [14]). But 
here we don't have a single maximum at x = ±1, but 
two light spots, i.e., two symmetric maxima around the 
positions of maximal downflow. They are caused by the 
concentration field. Its contribution to the shadowgraph 
intensity is opposite in sign to the temperature contribu- 
tion at the downflow roll boundaries for ip > 0, since 6 is 
negative [3] for all ethanol-water mixtures considered so 
far — for simplicity and better comparison we use here 
the same b as for the mixture with ip = —0.25. In SOC's 
with ip < 0, both contributions are positive, giving rise 
to the single maximum in Figs. 3d, 10b. 

Comparing a shadowgraph image in the Rayleigh 
regime (Fig. 12d) with one in the Soret regime (Fig. 12a), 
one sees that bright areas have become dark, and vice 
versa. This again reflects the facts that (i) concentra- 
tion and temperature fields contribute basically with op- 
posite signs to I(x,z), and that (ii) the former field is 
stronger than the latter field in the Soret regime, and 
vice versa high up in the Rayleigh regime. But the 
two bright spots for the smallest r-value in Fig. 12a at 
x = ±l,z « 0.75 indicate that concentration as well 
as temperature field contribute to the shadowgraph in- 
tensity, and that the concentration field is still strongly 
anharmonic. See also the lateral concentration profile in 
the fourth row of Fig. 12a. 

We should mention that in our calculations we did not 
wait for perfect stationarity of the SOC states. For the 
three states with the higher r, we observed a very slow 
relaxation process in the concentration field with a si- 
nusoidal spatial modulation given by the system length 
T = 20, while changes in the temperature and velocity 
field were nearly negligible. To filter out these long-time 
variations of SC, we presented in Fig. 12c an average of 
SC over 10 wavelengths. Further work is necessary to 
investigate whether this relaxation process reflects slow 
phase diffusion or weakly damped concentration modes. 



A reason for using localized traveling waves as initial 
conditions was to test if there exist TW's in the Rayleigh 
region for positive ip as predicted by the Galerkin model 
of Linz et al. [55] for free-slip impermeable boundary con- 
ditions. But in a more systematic search in the Rayleigh 
region for fluids with ip = 0.02 and ip = 0.05, we also did 
not find stable TW's. Furthermore, test simulations with 
free-slip impermeable boundary conditions did not show 
stable TW's for positive ip. On the other hand, Bigazzi 
et al. [22] found oscillatory behavior in their cylindrical 
convection cell which resembled TW's. Since oscillatory 
behavior in a cylindrical cell may have different causes, 
we think there is now more evidence against stable TW's 
for positive ip. Experiments in annular cells are most 
desirable to clarify this issue. 

We did not perform calculations right above the insta- 
bility despite the prospect of interesting wave number se- 
lection processes with possibly hysteretic character, since 
for our long system calculations are very time consum- 
ing, and since 3d effects have been predicted close to the 
instability [76,77,62]. 

VI. LEWIS- AND PRANDTL NUMBER 
DEPENDENCE 

In the preceding sections, we have discussed convective 
properties of binary mixtures as functions of the control 
parameters r and ip, while keeping the Lewis number L = 
0.01 and the Prandtl number a = 10 fixed. Yet L varies 
between 0.005 and 0.01 for ethanol-water mixtures in the 
temperature range between 10°C and 40°C [3], while a 
varies between 5 and 20 [3]. On the other hand, L = 0.03 
and a = 0.6 in 3 He- 4 He-mixtures, where experiments 
have also been performed. Furthermore, Steinberg et al. 
[30] report experiments with a benzene-methanol mixture 
at about L = 0.025, a = 7.5, ip = -0.045. 

We start the investigation of the L,cr, tp,r parameter 
space with the bifurcation diagrams of Fig. 13, each of 
which is for ip = —0.25. Characteristic properties for the 
four L — a combinations discussed in the following are 
given in Tab. II. Let us start from L = 0.01, a = 10 
(Fig. 2 and right curves of Fig. 13a, b) and increase the 
Lewis number to the 3 He- 4 He value L = 0.03, but keeping 
g — 10 constant (right curves of Fig. 13c, d). The Hopf 
frequency loh and stability boundary r osc [68] increase 
only slightly despite of the tripling of L. The saddle rf rw 
increases slightly, and the frequency at the saddle stays 
roughly constant. But r* decreases strongly, and thereby 
roughly halves the region r* — r^ w of stable TW's on the 
upper branch. But r* still lies slightly above r osc , so that 
the nonlinear TW at r osc has a small, but finite frequency 
w(r osc ) « JqU H - 

Let us now go back to our ethanol-water reference fluid 
(L = 0.01, a = 10) and reduce the Prandtl number to the 
3 He- 4 He value a = 0.6, while keeping L = 0.01 constant. 
For these parameters, r osc — 1 nearly halves, and wh de- 
creases to about | of the ethanol-water value (left curves 



of Fig. 13a, b). The saddle lies closer to r osc , so that 
the hysteretic TW region r osc —rf rw has become smaller. 
The frequency at the saddle has decreased to about \loh, 
compared to about \uh for the ethanol-water mixture. 
Since r* lies just above r osc , the region of stable TW's 
is significantly smaller than for ethanol-water with the 
higher Prandtl number. 

The two effects of (i) increasing L and (ii) reduc- 
ing a all in all "add up" for the 3 He- 4 He parameters 
L = 0.03, a = 0.6 (left curves in Fig. 13c, d). Note here 
that r* lies below r osc , i.e., the final state is nonlinear 
SOC. The region, r* — r^ w , of stable TW's on the up- 
per branch is smaller by a factor of 10 than for ethanol- 
water parameters, also if we scale with r osc — 1. But the 
qualitative structure of the bifurcation diagrams is the 
same for the four parameter combinations of Fig. 13. In 
particular there exists an upper TW branch. 

A more global picture of the changes of various bi- 
furcation properties with varying L and a is shown in 
the L — r plane of Fig. 14b and in the a — r plane of 
Fig. 15b. See also tables III and IV. With increasing L, 
the TW bifurcation threshold r osc first increases slightly 
and then strongly at L = C^IO" 1 ) (Fig. 14b). The Hopf 
frequency Wjj(fc = it), on the other hand, remains essen- 
tially constant for 0.001 < L < 0.1. The saddle r^ w lies 
roughly parallel to r osc in Fig. 14b, while r* decreases 
very strongly with increasing L. Thus, at L = 0.05, the 
transition to stationary convection is already below r osc 
(see also the frequency of the nonlinear state at r osc in 
Fig. 14a), and for L = 0.3 there is no longer a stable up- 
per TW branch, although r* has increased with increas- 
ing r osc . The vanishing of the upper TW branch here is 
similar to the vanishing in Fig. 9 for ip = —0.005. Note 
the rapid expansion of the region of stable TW's with 
decreasing L in Fig. 14b: the upper boundary r* shoots 
up from 1.65 at L = 0.01 to more than 5 at L = 0.005. 
This resembles the rapid increase of r* with decreasing 
tp around ip = —0.3 shown in Fig. 9b. The above de- 
scribed rapid variation with L presumably accounts for 
part of the differences that Liu and deBruyn [16] ob- 
served between their experiments at L = 0.005 and our 
calculations at L = 0.01. To understand the increase of 
the TW region with decreasing L, the boundary layer 
model of Bensimon et al. [51] might be a good starting 
point. However, in its present form, it fails to give the 
ending of the TW branch on the lower SOC branch and 
the extreme increase of r* , as we have seen in the L and 
%p variation. 

Let us now consider the variation with a while we keep 
L = 0.01 constant. With decreasing a, the Hopf fre- 
quency and r osc decrease, but below a « 0.14, the bifur- 
cation threshold r osc increases again. Also, the frequency 
of the nonlinear state at r osc decreases (Fig. 15a); indeed, 
it decreases faster than the Hopf frequency (see Tab. IV). 
In the considered parameter region, r* decreases down to 
the order of 1 with falling a like r* « a + &log(er). Down 
to a = 0.5, there exist still stable TW's above r osc , and 



stable TW's close to the saddle exist at least down to 
a = 0.1. Note that, in this small Prandtl-number regime, 
strong nonlinearities can also be seen in the velocity fields 
of TW and SOC states, which deserve further investiga- 
tion. 

Let us now address 3 He- 4 He mixtures with L = 0.03 
and a = 0.6 [36-45]. We have investigated their depen- 
dence on ip only in the range —0.25 < ip < —0.05. As 
we have discussed above, the qualitative structure of the 
bifurcation diagrams for ip = —0.25 is the same as for 
ethanol- water mixtures, but here r* lies just below r osc . 
This reduction of the region of stable TW's can also be 
seen in the variation with ip: while we still see stable 
TW's for ip = —0.1, there are no longer stable TW's for 
%p = —0.05. In comparison, in ethanol-water mixtures, 
there are stable TW's even for ip = —0.01. 

VII. CONCLUSION 

The main goal of this paper has been to provide a 
quantitative description of extended convection in bi- 
nary fluid mixtures and thus come to an understand- 
ing of the nonlinear stationary and traveling-wave states 
that occur as primary convective states near the convec- 
tive threshold. In particular, we wanted to elucidate the 
role of the concentration field. We should like to stress 
that the experimentally observable convective states on 
the upper stable parts of the bifurcation branches that 
we have investigated are strongly nonlinear. For exam- 
ple, for most of them, the concentration Peclet numbers 
wd/D = w/L measuring the magnitude of the nonlin- 
ear advective contribution to the concentration balance 
equation relative to the linear diffusive one were of order 
1000. The field structure of these states can be described 
quantitatively neither by a weakly-nonlinear expansion 
around their bifurcation threshold nor by simple, heuris- 
tic, one-mode amplitude-equation approaches. We there- 
fore solved the full hydrodynamic field equations numer- 
ically with realistic boundary conditions, using an ex- 
plicit finite-differences method. We restricted ourselves 
to 2d convection in the form of straight rolls by suppress- 
ing spatial variations along the roll axes. In most cases, 
we enforced the wavelength of the states, by imposing 
laterally periodic boundary conditions, to be A = 2, i.e., 
close to the wavelengths seen in experiments. 

Typical bifurcation behaviour for negative Soret- 
coupling — In a bifurcation diagram of convective inten- 
sity versus reduced Rayleigh number, an unstable TW so- 
lution branch bifurcates backwards at a Hopf bifurcation 
threshold r osc out of the conductive state. This TW 
branch turns over at a saddle r s TW below r osc and be- 
comes stable. The upper stable branch ends by merging 
at r* above r^ w into the SOC solution branch. Moving 
along the TW branch, the TW frequency continuously 
decreases from the Hopf frequency at threshold down to 
zero at the end point r* . The SOC branch bifurcates at 



Tstat backwards out of the conductive state or is discon- 
nected from it after r sta t has reached infinity. This lower 
branch bends forwards at a saddle r s soc . By suppress- 
ing TW convection with phase pinning boundary condi- 
tions, we determined the upper SOC branch all the way 
down to its saddle. However, allowing for phase propa- 
gation, SOC's are unstable below r*: for r < r^ w they 
start propagating and decay into the conductive state, 
and in the interval r s TW < r < r* , they start propagat- 
ing and end up as stable TW's. Viewed from the SOC 
branch above r* , this branch loses stability in a mirror- 
symmetry-breaking bifurcation to a stable left or right 
TW whose frequency grows ~ (r* — r)s close to the 
bifurcation threshold r* , in good agreement with experi- 
ments in ethanol-water mixtures in an annular cell [5,11]. 
These annular cells are closed convection channels with 
a rectangular cross-section. Apart for the channel walls, 
there are no lateral walls which hinder the traveling of 
the TW's as in laterally finite channels. The transition 
TW <-> SOC at r* is not forced by lateral walls but is 
an intrinsic effect which was described by a concentra- 
tion boundary layer expansion [51] around the states of 
steady convection in a pure fluid. This model of Bensi- 
mon et al. [51] can describe this transition and shows the 
connection of the upper TW branch via the lower TW 
branch with the conductive state at r osc . 

Structure of the fields — Both the velocity and temper- 
ature fields of SOC's and of TW's look similar in struc- 
ture and in strength to those in a comparable stationary 
state in pure-fluid convection. But in the TW state of the 
mixture, the temperature wave is slightly phase-shifted 
relative to the velocity wave. The shift varies nearly lin- 
early with the TW frequency as predicted earlier [72,55]. 
The lines of constant phase of the first lateral TW Fourier 
modes being curved, the waves of temperature, concen- 
tration and velocity are not planar. 

While velocity and temperature waves are nearly har- 
monic, the concentration wave has a strongly anhar- 
monic, trapezoidal profile. The latter shows that, in a 
TW, the concentration is almost homogeneous within 
each roll, but at alternatingly high and low concentra- 
tion levels, with a linear variation between adjacent rolls. 
This structure can best be understood within the refer- 
ence frame that comoves with the TW's phase velocity 
v p , so that all fields are stationary therein. With the dif- 
fusion constant being small, concentration is transported 
predominantly passively by the velocity field. Thus, one 
can qualitatively infer from its streamlines in the co- 
moving frame the concentration distribution in a TW. 
Because of the contribution — v p z to the stream func- 
tion from the phase velocity, the roll-like regions with 
closed streamlines are shifted alternatingly towards the 
top and bottom plates, and in addition there are open 
streamlines. They meander between the roll-like regions 
and separate them from each other and from the plates, 
in such a way that the boundary layer with the Soret- 
induced concentration surplus (deficiency) at the top 
(bottom) plate feeds, for instance, in a right TW only 



the right (left) turning roll regions: high (low) concen- 
tration plumes emanating from the top (bottom) bound- 
ary layer are bent into the right (left) turning roll re- 
gions where the fluid is mixed up. So the concentration 
wave profile resulting from this combination of bound- 
ary layer behaviour perpendicular to the open velocity 
streamlines and mixing behaviour within the regions of 
closed streamlines cannot be described by a few lateral 
Fourier modes, e~ mkx , as the velocity or temperature 
wave. This clearly shows that stable TW's on the upper 
bifurcation branch are strongly nonlinear states whose 
structural dynamics differs significantly from marginally- 
stable linear TW modes at the Hopf bifurcation thresh- 
old. 

The phase velocity or frequency of the TW is the sin- 
gle most relevant parameter that governs the structural 
properties of TW's. With decreasing v p , the streamlines 
of the velocity field in the comoving frame become more 
and more mirror symmetric around up- and downflow 
positions, the regions of closed (open) streamlines grow 
(shrink) , their alternating up- and downwards shift from 
the middle of the layer is reduced, the feeding of concen- 
tration into the different roll regions becomes more sym- 
metric, and the concentration contrast between adjacent 
rolls — and with it the plateau-height in the trapezoidal 
profile of the concentration wave — is reduced until at 
r* the SOC concentration distribution is perfectly mirror 
symmetric around up- and downflow. Then concentra- 
tion plumes ascending (descending) at the positions of 
maximal up (down) flow from the bottom (top) bound- 
ary layer inject low (high) concentration symmetrically 
into left- and right-turning rolls, where the fluid is well 
mixed in the bulk of the layer. The symmetric mixing 
operation of SOC's homogenizes the concentration more 
efficiently than the asymmetric mixing at different lev- 
els in adjacent TW rolls. This can be quantified by our 
mixing parameter M measuring the rms variations of the 
concentration field. 

The plume-like transport of concentration out of the 
boundary layers near the plates into the bulk of TW's 
and SOC's is a consequence of the smallness of the dif- 
fusion constant and thus is similar to thermal plumes in 
high Rayleigh number convection, where convective heat 
transport dominates in the bulk over diffusive transport. 

Symmetry properties — TW's and SOC's show a mir- 
ror glide symmetry under lateral translation by half a 
wave-length combined with reflection through the hori- 
zontal midplane of the layer. This can easily be seen, e.g., 
in side-view shadowgraphs [15] which, by the way, yield 
detailed information on the concentration distribution in 
the fluid. 

Mean transport of TW's — TW's cause a lateral mean 
flow U(z) that is several orders of magnitude smaller than 
v p and that is related [55,54] to the ^-variation of the 
phase of the first lateral Fourier mode of the lateral ve- 
locity field. On the other hand, and more importantly, 
the phase shift between vertical velocity and tempera- 
ture (concentration) wave generates a significant mean 



lateral convective heat current (u ST) (concentration cur- 
rent (uSC)) that flows parallel (opposite) to the TW 
propagation in the upper half of the fluid layer and vice 
versa in the lower half, so that the net currents vanish. 
Both current amplitudes roughly grow proportional to 
v p . The lateral convective transport of heat and concen- 
tration by these currents is different from and unrelated 
to the mean flow U(z). 

Changing the Soret coupling — The upper (lower) exis- 
tence limit r* (r^ w ) of TW's increases strongly (slightly) 
as the coupling parameter %p becomes more negative. 
On the other hand, as the magnitude of the negative 
separation ratio ip is reduced, the TW end point r* 
moves downwards on the SOC branch, so that r* ap- 
proaches the saddles r s soc and r^ w . The upper part of 
the TW branch shrinks to zero when r* and r^ w meet. 
That happens close to r s soc . In ethanol- water mixtures 
(L = 0.01,cr = 10), we found for -ip = -0.01 still a small 
r-range of stable TW's, cf. [49, Fig. 5.4], but not for 
tP = -0.005. So we conclude that for tp > -0.005 the 
TW endpoint has already moved around the SOC sad- 
dle and lies on the lower SOC branch, while the upper 
SOC branch has become stable all the way down to r s soc . 
A negative Soret coupling that is smaller in magnitude 
implies smaller concentration variations with smoother 
profiles, broader boundary layers, and wider plumes. 

Positive Soret coupling — Calculations done for ip = 
0.1 in a system of lateral periodicity length T = 20 
showed in the Soret regime, r < 1, but well above the in- 
stability and also in the Rayleigh regime, r > 1, always an 
evolution into SOC states with wavelength A = 2. This 
should be compared with experiments in narrow convec- 
tion channels that suppress squares and favour rolls. We 
used localized traveling waves as initial conditions to test 
if there exist stable TW's at -0 = 0.1 but found none 
— also not with free-slip impermeable boundary con- 
ditions. The characteristic structure of the SOC con- 
centration field — boundary layers, plumes, and mixing 
within the rolls — appears as for negative Soret coupling. 
Here, however, there is a concentration surplus at the 
lower, warm plate. Reducing r into the Soret regime, the 
flow intensity decreases, diffusive transport contributes 
more and more, concentration profiles become smooth, 
and boundary layers and plume structures become broad. 
But SOC's in the Soret regime are still nonlinear states 
as documented by the mixing parameter M being well 
below the conductive reference value M con d = 1. 

Lewis- and Prandtl number dependence — Our calcu- 
lations at ip = -0.25 done for 0.005 < L < 0.3 and 
0.1 < a < 10 showed that, upon decreasing L or increas- 
ing a beyond the reference values L = 0.01, a = 10 of 
ethanol water, the existence range of stable TW's on the 
upper branch widens with increasing r*. On the other 
hand, increasing L or decreasing a eventually pinches off 
their existence range as r* and rf rw come together. But 
otherwise the bifurcation structure of TW's and SOC's 
remains unchanged for negative Soret coupling. 
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TABLE II. Linear and nonlinear properties of the bifur- 
cation diagrams in Fig. 13 for tp = —0.25 are given for four 
combinations of L and a. 
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TABLE IV. <T-dependence in mixtu- 

res with L = 0.01, i\> = —0.25. Given are r osc (k = tt) and 
cuh(& = 7r) [68]. The last four columns contain properties 
of the nonlinear convective states on the upper bifurcation 
branches at r osc . 
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TABLE III. Independence in mix- 

tures with g = 10, tp = —0.25. Given are r osc {k = tt) and 
Wif(t = n) [68]. The last four columns contain properties 
of the nonlinear convective states on the upper bifurcation 
branches at r osc . 
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FIG. 1. ^-dependence of linear stability properties of the 
heat conducting state of ethanol- water mixtures [68]. (a) Crit- 
ical reduced Rayleigh numbers r osc (dashed line) of the oscil- 
latory and r s tat (full and dash-dotted line) of the stationary 
instabilities, (b) Critical reduced wave numbers of the oscil- 
latory (dashed line) and stationary (full line) instability for 
r > 0. The critical wave number of the stationary instability 
at r < vanishes, (c) Critical Hopf frequency ujh- 



FIG. 2. Bifurcation diagrams of extended convective states 
with wavelength A = 2. (a) Frequency, (b) mixing parame- 
ter M, (c) Nusselt number. The TW solution branch bi- 
furcates backwards (this unstable branch is not shown) at 
r OS c = 1.3347 with frequency u>h = 11.235 [68], becomes sta- 
ble (full circles) at r^w = 1.213, and ends at r* = 1.65. SOC's 
(open squares) are stable only above r* . By phase pinning, 
they can be stabilized (open triangles) down to r s soc = 1.08. 
The lower unstable SOC solution branch (not shown) is dis- 
connected from the conductive solution and lies almost par- 
allel to the r-axis. The dotted line in (c) shows the Nusselt 
number in a pure fluid with a = 10. 



FIG. 3. Snapshots of three representative TW's traveling 
to the right and a SOC state for L = 0.01, a = 10, V = -0.25. 
(a) TW at r = 1.2155 close to the saddle with uj « 0.3loh, (b) 
TW at r = 1.2459 with uj « \uj h , (c) TW at r = 1.4181 with 
uj « 0.07u;h, and (d) SOC state at r = 1.8232. See Fig. 2 
to locate the r-values in the bifurcation diagrams. The wave- 
length is A = 2, and x = denotes the position of maximal 
upflow. The colour coded temperature field ST (first row), 
with arrows denoting the velocity field u, is shown only for 
the fastest TW (a) and the SOC state (d). The second and 
seventh row show the sideview shadowgraph intensity I(x, z) 
(3.15) with b = —0.919 in grey scales. The concentration field 
SC in the third and sixth row is given in the same colour code 
for all four states. Red, green, and blue represent low, mean, 
and high concentration, respectively. The structure of the 
concentration field can be related to the dotted streamlines 
of the velocity field u in the frame comoving with the phase 
velocity of the states. Fourth and fifth rows show lateral pro- 
files of w (thin line), 40 ST (triangles), and 400 SC (squares) 
at midheight, z = 0.5. 



FIG. 4. Vertical structure of TW and SOC states. Columns 
show from left to right the states of Fig. 3: (1) TW at 
r = 1.2155, (2) TW at r = 1.2459, (3) TW at r = 1.4181, 
(4) SOC state at r = 1.8232. Row (a) shows ^-dependence 
of lateral Fourier modes 8o and \6 n \ (n = 1, 2, 3) of the tem- 
perature field. In row (b), the corresponding modes of the 
concentration field, SC, are presented. For better resolution 
especially at the plates, Fourier modes were evaluated also 
vertically halfway between grid positions. Row (c) shows ver- 
tical profiles of the concentration field SC at ^-positions A 
= 0.35, B = 0.85 in the right turning roll and A' = —0.65, B' 
= —0.15 in the left turning roll differing by -|. Maximal up- 
flow is at x = as in Fig. 3. For comparison, the conductive 
profile is shown as a dotted line. 



FIG. 5. Rayleigh-number dependence of lateral Fourier 
modes \ f n \ of: (a) vertical velocity / = w, (b) temperature 
/ = 8, (c) concentration / = SC for TW (closed symbols) 
and SOC states (open symbols). SOC's below r* are stabi- 
lized by phase-pinning boundary conditions. Presented are 
at z = 0.5 (dots for TWs and squares for SOC's), |/ 2 | at 
z = 0.25 (lozenges), and | jTs | at z = 0.5 (triangles). For com- 
parison, the corresponding modes of SOC states in a pure 
fluid, ip = 0, are included in (a) and (b) by lines: full for 
n = 1, dotted for n = 2, and dashed for n = 3. 



FIG. 6. Frequency dependence of lateral Fourier modes 
of the concentration field of TW's at different z positions, 
(a) \SCi\, (b) phase shift between concentration wave and 
wave of vertical velocity, (c) lateral average SCo- In (c), 
SCo(z = 0.5) = is not included. The highest frequency 
denotes the TW nearest to the saddle, r?p W , while u> = at 
r*. For u> — > 0, the phase shift in (b) approaches -k in the 
bulk (z = 0.5 and z = 0.25) but drops to zero at the plates 
(z = 0). 



FIG. 7. Streamlines in the x-z plane for time-averaged 
transport in the right-going TW of Fig. 3a. (a) mean flow 
(u), (b) heat current (Q) — Q CO nd, and (c) concentration cur- 
rent (J). 



FIG. 8. Lateral transport of TW's. (a) Amplitude of lat- 
eral flow \U(z = 1/2) | . (b) Maxima of the mean lateral cur- 
rents of heat, (Qx) (triangles), and of concentration, (J x ) (cir- 
cles), (c) Shows the TW frequency uj for comparison. 



FIG. 9. (a) Frequency of nonlinear convective states on 
the upper bifurcation branch at r oac reduced by the Hopf fre- 
quency ujh ■ TW's are shown by full dots and SOC states with 
to = as full squares. See also Tab. I. (b) Double-logarithmic 
phase diagram. Symbols come from simulations done with 
prescribed wavelength, A = 2. The conductive state is sta- 
ble below and unstable above the TW bifurcation threshold 
r sc(k = 7r) (thin line) [68]. Stable TW's on the upper branch 
exist in the shaded area between the TW saddle r^ w (full 
dots) and r* (full squares), where the TW branch merges into 
the SOC branch. SOC states are stable above r* or — for i\> 
values close to zero where we have not found stable TW's — 
above the SOC saddle r%oc (line with filled squares in the 
lower right corner). Where stable TW states are found, SOC 
states are unstable below r* . But the upper SOC branch can 
be stabilized down to the SOC saddle r% oc (open squares) by 
suppressing TW's. At i\> = —0.01, where r s soc « r^w ~ r* , 
there still exists a very small area of stable TW's, while for 
ip = —0.005 we have not found any stable TW. 



FIG. 14. (a) Frequency of nonlinear states at r osc on the 
upper bifurcation branch versus Lewis number. Note the log- 
arithmic scale of L. (b) L-r phase diagram. Symbols denote 
final states obtained in simulations. The dashed curves r* and 
r^w denoting the borders of the dashed area of stable TW's 
are hand-drawn interpolations. The thin curve, r oac [68], is 
the TW bifurcation threshold. Some properties of the phase 
diagram are given in Tab. III. 



FIG. 15. (a) Frequency of nonlinear states on the upper 
bifurcation branch at r osc versus Prandtl number. Note the 
logarithmic scale of a. (b) a — r phase diagram. Symbols 
represent final states obtained in simulations. The dashed 
curves, r* and r s TW , denoting the borders of the dashed area of 
stable TW's are hand-drawn interpolations. The thin curve, 
r sc [68], is the TW bifurcation threshold. Some properties 
of the phase diagram are given in Tab. IV. 



FIG. 10. Snapshots of convective states 
at ?/> = -0.01, L = 0.01, a = 10. (a) TW at r = 1.01483 
with u « 0.074, (b) SOC state at r = 1.0212. The presen- 
tation is as in Fig. 3. However, to resolve the concentration 
variation, its colour code is different as indicated. The lateral 
profiles at z = 0.5 in the last row show w (thin line), 20 ST 
(triangles), and 3000 SC (squares). 



FIG. 11. Bifurcation diagrams for positive Soret coupling. 
The simulations (squares) of a system of length T = 20 pro- 
duced stable SOC states. Their wavelength was always A = 2. 
(a) Mixing parameter M, (b) Nusselt number. The SOC 
branch in a pure fluid is shown as the dotted curve. Thresh- 
olds are r s tat(k c = 0) = 0.04216 and r s tat(k = tt) = 0.06012 
[68]. 



FIG. 12. Snapshots of the first four SOC states of Fig. 11. 
The states at (a) r = 0.506 and (b) r = 0.9116 lie in the Soret 
region. Those at (c) r = 1.1142 and (d) r = 2.026 lie in the 
Rayleigh region. The lateral profiles at z = 0.5 show w (thin 
line), 40 ST (triangles), and 1000 5C (squares). Otherwise the 
presentation is as in Fig. 3. 



FIG. 13. Bifurcation diagrams of frequency and Nusselt 
number for ip = —0.25 and four combinations of L and a over 
a common interval of Rayleigh numbers. See also Tab. II. 
Stable TW's are presented by dots and SOC states by squares. 
Curves are guides to the eye. Unstable TW branches (dashed 
lines) connecting to ujh in the upper row and to r osc (crosses) 
in the lower row are schematic. SOC states in the pure fluid, 
ip = 0, are shown by dotted lines in the lower row. 
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